X = 0:1000:50000;
Y = 0:1000:70000;
aa = length(X);
b = length(Y);
t= 0:1:360;
map =zeros(aa,b,2);
L1=4000;
plan=zeros(1,7*5);
M=70000;
%% 开始循环过程
for a= 0.3:0.1:0.9
for i=1:aa
for j=1:b
    p=0;
    for V1=70:12.5:95
        for V2=80:10:100 
             [t,xxyy] = ode45( @(t,y)fly_the_air(t,y,L1,V1,V2,a,M),t,[X(i),Y(j),50000,50000,50000,20000]);
            [status ,tt] = check_status(xxyy,M);

%  status 的状态
%  0 - 未知
% 1 - 拦截成功
% -1  ： 曲率不满足条件
% 2 ： 突围成功
            if status ==1 
                p=1;
            end
            if status==2
                p=0;
                break;
            end
        end
        if p==0
     disp(i)
     disp(j)
     disp("----")
            break
        end
    end
    if(p==1)
        map(i,j)=1;
    end
end
end
end

            